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ANALYSIS AND COMPENSATION OF AN AIRCRAFT SIMULATOR 
CONTROL LOADING SYSTEM WITH COMPLIANT LINKAGE 

By Paul R. Johansen* and Richard E. Bardusch 
Langley Research Center 

SUMMARY 

A hydraulic control loading system for aircraft simulation was analyzed to find the 
causes of undesirable low frequency oscillations and loading effects in the output. The 
hypothesis of mechanical compliance in the control linkage was substantiated by compar- 
ing the behavior of a mathematical model of the system with previously obtained experi- 
mental data. 

A compensation scheme based on the minimum integral of the squared difference 
between desired and actual output was shown to be effective in reducing the undesirable 
output effects. The structure of the proposed compensation was computed by use of a 
dynamic programing algorithm and a linear state space model of the fixed elements in 
the system. 


INTRODUCTION 

Control systems in aircraft which convert stick, wheel, and pedal motions into con- 
trol surface deflections constitute the interface between pilot and airframe. In addition 
to the pilot input function, such systems provide force cues to the pilot which assist him 
in assessing the extent and rate of controlled maneuvers. The aspect of an aircraft con- 
trol system which concerns these cues is called control loading. 

One objective of an aircraft flight simulator is to familiarize pilots with the control 
forces required to perform well-executed maneuvers. Thus, the control loading system 
in an aircraft simulator must be capable of accurately reproducing force cues typical of 
the real aircraft. Most simulator control loading systems in current use consist of an 
electrohydraulic position servo driven by an analog computer, which can generally be 
programed to simulate a wide spectrum of aircraft characteristics. 

The requirement of fast response over a wide performance range, together with 
inertial and elastic effects in the control linkage, can cause unsatisfactory performance of 
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a control loading system in some ranges of operation. The purpose of this study was to 
examine a particular wide-range control loading system in order to determine the causes 
of unsatisfactory performance and to devise an alternative control scheme which would 
improve the behavior of the system over the specified range of operation. 

The problem was attacked in two phases. The analysis phase of the study was 
needed to determine the dynamic characteristics of the fixed components of the system, 
that is, the plant. Prior experimental data showed that there was an undesirable low 
frequency oscillation in the elevator system which could not be attributed to known ele- 
ments of the system. The hypothesis of mechanical compliance in the control linkage 
was substantiated by comparing the behavior of a mathematical model of the system with 
previously obtained experimental data. 

The design of an alternative control scheme commenced when the mathematical 
model of the fixed plant, including linkage compliance, was judged to be capable of sat- 
isfactorily predicting the behavior of the real plant. The form of the proposed compen- 
sation was computed by use of a dynamic programing algorithm and a linear state space 
model of the plant. Input and feedback control coefficients were obtained for the entire 
simulation range of the system by minimizing a quadratic cost index. Behavior of the 
closed-loop system, as predicted by the revised mathematical model, indicated that the 
low frequency structural oscillation can be reduced. The compensated system was shown 
to be relatively insensitive to variations in the control coefficients, and certain feedback 
paths could be eliminated without seriously degrading the performance of the system. 

The revised model indicates that the control scheme is physically feasible. 

SYMBOLS 

Values are given in both SI Units and U.S. Customary Units; the measurements and 
calculations were made in U.S. Customary Units. 

A cross-sectional area of actuator piston, cm2 (in2) 

B modeled damping coefficient, N-cm-sec (lb-in-sec) 

c damping coefficient of linkage, N-sec/cm (lb-sec/in.) 

c control coefficient vector 

c^jCg,. . .,c ? components of control coefficient vector 

Cj valve current gain, cm 3 /sec-mA (in 3 /sec-mA) 
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c 


p 



I 

J 

k 


K 2 

K 3 


P 


valve pressure gain, cm5/N-sec (inVlb-sec) 

plant-input control variable, V 

output of load cell circuit, V 

output of potentiometer circuit, V 

input test signal, V 

load cell force, N (lb) 

pilot input force, N (lb) 

modeled spring rate, N-cm/rad (lb-in/rad) 

maximum value of modeled spring rate, N-cm/rad (lb-in/rad) 

current, mA 

rated valve current, mA 

modeled control inertia, N-cm-sec^ (lb-in-sec2) 

mass moment of inertia of column, N-cm-sec 2 (lb-in-sec 2 ) 

spring rate of linkage, N/cm (lb/in.) 

transducer gain, V/N (V/lb) 

admittance of valve circuit, mA/V 

gain of potentiometer circuit, V/cm (V/in.) 

distance between link attachments on bellcrank, cm (in.) 

distance between piston link attachment and bellcrank pivot, cm (in.) 

load pressure, N/cm 2 (lb/in2) 
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Pi 

P2 

P r 

Ps 

q 


<*r 

R 1’ R 2 

Ru 


pressure on side 1 of piston, N/cm 2 (lb/in 2 ) 
pressure on side 2 of piston, N/cm 2 (lb/in 2 ) 
valve return pressure, N/cm 2 (lb/in 2 ) 
valve supply pressure, N/cm 2 (lb/in 2 ) 
fluid flow rate, cm 3 /sec (in 3 /sec) 

component of linearized valve flow rate, cm 3 /sec (in 3 /see) 
rated valve flow rate, cm 3 /sec (in 3 /sec) 


weighting factors 


lever arm of column base, cm (in.) 


lever arm of pilot force, cm (in.) 


t 

T 

u 


time, sec 


column torque, N-cm (lb-in.) 


input vector 


performance index 


E 


effective volume of compressed fluid, cm 3 (in 3 ) 


state vector 


displacement of column base, cm (in.) 
displacement of actuator piston, cm (in.) 


y output vector 

[a] ,[b] , [c] ,[D] , [e] matrices 
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integer 
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DEVELOPMENT OF MATHEMATICAL MODEL 


A complete aircraft simulator control loading system includes systems which pro- 
vide force cues in the elevator, aileron, and rudder control modes; however, only an ele- 
vator system is considered here. Ordinarily, although elevator, aileron, and rudder sys- 
tems are independent of one another, they are similarly structured and, therefore, an 
analysis of the elevator system would be applicable to aileron and rudder systems with 
slight modifications. In addition, the pilot control elements associated with the elevator 
system vary from one aircraft to another more widely than do the similar elements for 
the aileron and rudder systems; consequently, the elevator system would represent the 
loading mode most likely to present difficulties in the design of a wide-range simulation 
system. 

The elevator control loading system examined here is shown in mechanical sche- 
matic form in figure 1. The control loading system senses the pilot input force by means 
of a load cell in the actuator-column linkage and transmits an electrical signal e F to 
the computer which in turn commands the positioning system. The difference between 
the computer output and the position feedback ep controls the fluid flow and thereby the 
motion of the control column. The computer can be programed to provide a wide range 
of dynamic characteristics typical of control loads which occur in real aircraft; however, 
programable nonlinear effects, such as breakout, hysteresis, and velocity limiting, are 
not considered herein. 

As figure 1 indicates, the control loading system comprises a doubly closed-loop 
system whose input force and output position are transmitted through the same element, 
that is, the pilot's control column. Ideally the inner loop, which contains the servovalve 
and actuator, should have a fairly flat frequency response over a wide bandwidth. Rapid 
response is desired, but a change in phase margin must be made to provide stability when 
the outer loop is closed. Thus, the inner loop (hydraulic servo) is designed to be some- 
what overdamped. A development of the mathematical model of the original system 
follows. 


Hydraulic Servo Model 

The principal components of the hydraulic servo are the amplifier, valve, actuator, 
load, and position transducer. The amplifier is assumed to have sufficient bandwidth to 
be considered a constant multiplier whose input is a voltage signal and whose output is a 
valve current. The servovalve is a four -way, two- stage electrohydraulic control valve 
whose input is a current signal and whose output is a fluid flow rate through the actuator. 
The static flow-pressure characteristics of the valve are essentially parabolic, typical of 
orifice flow-pressure relations (ref. 1). This relationship 
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( 1 ) 



which is illustrated in figure 2, may be linearized about a nominal operating point through 
the use of a modified Taylor's series (ref. 2) as follows: 


q = c ii - C p p 

where 


( 2 ) 


Ci 


3q 


i=i 

P=P 


C P = 


Oq 

Op 


1=1 

P=P* 


and i* and p* are the valve current and load pressure at the nominal operating point 
In this model, i* was chosen to be 20 percent of rated current and p* was chosen to 
be 33 percent of supply pressure. Servovalve frequency response (that is, the ratio of 
load flow with zero load pressure drop to input current) was assumed to be similar to a 
linear second-order system with a natural frequency of 60 Hz and damping ratio of 0.9 
(ref. 3). The resulting equations describing the dynamic behavior of the servovalve are 


\ + + w n 2q i = C i w n 2i 

and 

q = q A - C p p (4) 

Thus, the flow to the actuator is expressed as a sum of linear terms in control current 
and load pressure. 

The dynamic equation governing the motion of the actuator can be determined by 
considering the flow continuity equation (ref. 1) or 


q = Ax p + (V E /4^p 


( 5 ) 



Although leakage between the piston and cylinder is often included in continuity equations 
for actuators, it is neglected in this analysis because the actuator piston seal in this sys- 
tem reduces leakage flow to an extremely small fraction of displacement and compress- 
ibility flows. 


Load and Linkage Dynamics 

A peculiar characteristic of a simulator control loading system is that the input and 
output elements are one and the same. The pilot input to the system is a force applied to 
the control column, and the output is the angular displacement of the column. A sche- 
matic diagram of this input-output arrangement is shown in figure 3. Two models will be 
formulated: (1) a rigid model based upon the assumption that all elements of the linkage 

and supporting structures are rigid, and (2) a compliant model based upon the assumption 
that the input-output linkage contains flexible elements which may be considered as one 
member in the analysis. 

Rigid linkage model .- The equation for the motion of the actuator piston and control 
column may be formulated by applying Newton's second law to a free-body representation 
of the control column. The torques about the pivot point of the column (see fig. 3) are 
summed to give 

/_j T — RpFp + aRjjAp = JO (6) 


From the geometry of the column, 0 = x^/r^ (for small motions), where x^ is the 
linear displacement of the base of the column. Thus, equation (6) may be written as 

x b = (aR b 2 A/j)p + (R b Rp/j)F p (7) 

Expressed in terms of the actuator displacement Xp, equation (7) becomes 



The mass of the actuator piston and load cell is assumed to be negligible compared to the 
inertia of the column. 

Compliant linkage model .- The dynamic equations which describe the motion of the 
actuator piston and control column can be determined from a free-body representation of 
the load cell attached to the actuator output. In this case, the linkage is assumed to con- 
tain flexible elements, and the dynamic properties may be contained in the linkage mem- 
ber immediately connected to the load cell (see fig. 4); that is, 
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( 9 ) 


£ F C = A P = k ( x p - ax b) + c ( x p - ax b) 

Equation (9) may be arranged in the following manner: 

x p = - (k/ c)x p + (A/c)p + (a!k/c)x b + ax b (10) 

Since the load cell and piston are assumed to contribute no inertial force, the equation of 
motion for the column is given as before by equation (7). The spring rate and damping 
coefficient can be computed from experimental observations of column motion and know- 
ledge of the column inertia. 

It is assumed that the displacement of the actuator piston and column can be mea- 
sured directly with transducers whose dynamic contributions are negligible and also that 
the load cell contributes no dynamic effects. 

Analog Computer 

The analog computer receives force signals from the load cell and generates dis- 
placement voltages which drive the hydraulic servo. In order to reproduce the loading 
effects typical of a wide variety of aircraft controls, the computer may be programed to 
simulate a linear second-order system with variable natural frequency and damping ratio. 

The dynamic equation solved by the analog computer is 

I6 + B6+G6 = T (11) 

The input and output of the analog computer are scaled voltages, proportional to the actual 
column torque and desired output, respectively. Accordingly, equation (11) may be 
written as 


16 + B6 + G6 = G max K 1 Ap (12) 

The hydraulic servo input is related to the analog computer output voltage by 

i = K 2 (6-K 3 x p j (13) 

The effect of compliance in the load linkage can now be assessed. The presence of 
a compliant member introduces an additional degree of freedom and associated phase 
lag; consequently, the rigid linkage assumption leads to incorrect predictions of system 
performance and inappropriate feedback policies. Inclusion of a compliant linkage mem- 
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ber in the mathematical model of the system enables the model to predict response phe- 
nomena similar to that observed in the laboratory. Furthermore, the hypothesis of 
unavoidable compliance in the control linkage provides a basis for compensation of the 
system. 


State Vector Models 

Linear state vector formulation is a convenient means by which analysis, simula- 
tion, and compensation of the control loading system may be accomplished. A typical 
form of these vector equations is 


y-M 


X + 


X + 


[b]G 

[d]G 


(14a) 


(14b) 


where [a], [b], [c], and [d] are constant matrices. 


The complete dynamic model of the control loading system with rigid control link- 
age is shown as a block diagram in figure 5. The numerical values used in this system 
are given in the appendix. The appropriate assembly of equations (3), (4), (5), (8), (12), 
and (13) which expresses the rigid model in the form of equations (14a) and (14b) becomes 



Xp 


0 

1 

0 

0 

0 

0 

0 

X P 


0 


^p 


0 

0 

a 2 R b 2 A/j 

0 

0 

0 

0 

x p 


aR^Rp^J 


p 


0 

-4A/yV E 

-^c p /v E 

4£/v e 

0 

0 

0 

P 


0 

d 

dt 

% 

“ 

0 

0 

0 

0 

1 

0 

0 


+ 

0 



'Ci„ n 2 K 2 K 3 

0 

0 

V 


Ci^nS 

0 

% 


0 


6 


0 

0 

0 

0 

0 

0 

1 

6 


0 

i 

6 


0 

0 


0 

0 

-G/I 

-B/I 

J 

6 


0 


9 = 


<xR b 


F p (15) 


(16) 


Figure 6 shows a block diagram of the model having compliant linkage. The 
numerical values used in this system are given in the appendix. The complete dynamic 
model including compliant linkage is written in state vector form by using equations (3), 
(4), (5), (7), (10), (12), and (13); that is, 
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* - 
X P 


-k/c 

A/c 

0 

0 

a 

0 

ark/c 

a 1 

x p 


- 0 




4Aj3k/v E c 

-4^(a 2 + C p c)/v E c 

40/v e 

0 

0 

0 

-4A0ak/V E c 

-4A0a/v E 

p 


0 


li 


0 

0 

0 

1 

0 

0 

0 

0 



0 

d 



“ C i w n 2K 2 K 3 

0 

-<0 2 
n 


C i“n 2K 2 

0 

0 

0 


+ 

0 

dt i 

5 


0 

0 

0 

0 

0 

1 

0 

0 ! 

a 


0 


6 


0 


0 

0 

-G/I 

-B/I 

0 

0 : 

6 


0 


*b 


0 

0 

0 

0 

0 

0 

0 

1 

x b 


0 


x b 


0 

aR b 2 A/j 

0 

0 

0 

0 

0 

0 

k b 


R b R p /j 


( 17 ) 



(18) 


The input signal e T , shown at the input to the valve amplifier in figures 5 and 6, 
represents a voltage signal used for obtaining system frequency response data both in the 
laboratory (see fig. 7) and in the mathematical models. 


COMPARISON OF RIGID AND COMPLIANT MODELS 


All numerical values for coefficients in the rigid linkage model were determined 
from component manufacturers' specifications, geometric considerations, and analytical 
calculations. Simulation of the system using this model revealed that such a model is 
incapable of predicting the low frequency resonances observed in the laboratory. 

Inclusion of a compliant linkage member between the actuator piston and the control 
column in the mathematical model provided the means by which approximate laboratory 
behavior could be predicted. The linkage spring rate and damping coefficient were chosen 
so that the observed resonance in the hydraulic servo was matched by the mathematical 
model. Predicted open-loop and closed- loop servo frequency responses are shown in fig- 
ures 8 and 9, respectively. Resonance peaks, which are seen to occur at approximately 
12 Hz with a damping ratio of 0.1, correspond favorably with laboratory observations. 

The overall system open-loop frequency response, shown in figure 10, illustrates 
the sizable decrease in phase margin introduced by compliance in the load linkage. This 
difference in phase margin is a plausible explanation for the discrepancy between behavior 
predicted by the rigid model and that observed in the laboratory; the closed-loop system 
becomes unstable for low values of computer- modeled spring rate G which correspond 
to high loop gain. Typical predicted closed-loop responses to step pilot forces for the 
rigid and compliant models can be compared in figures 11(a) and 11(b). The response of 
the compliant model displays the characteristics of the laboratory system and thus quali- 
fies it for use in developing compensation for the laboratory system. 
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COMPENSATION OF THE SYSTEM 


The problems of linkage compliance and resulting output oscillations have at least 
two possible solutions. One solution is the removal of sources of compliance in the link- 
age by altering or stiffening the linkage members, linkage bearings, and bearing mounts. 
This method, however, may be unsatisfactory because of additional weight or space 
requirements, and quite likely such a solution would, at best, elevate the structural 
resonance rather than eliminate its effects. 


Another solution, which is developed here, considers the linkage compliance to be 
an inextricable component of the fixed plant. Accordingly, the compensation proposed in 
this section requires knowledge of the resonant frequency of the linkage. This informa- 
tion is obtained from the laboratory system by measuring the piston and column responses 
in tests of the loaded hydraulic servo. 


The fixed plant, for purposes of compensation, is shown in figure 12. The plant is 
considered to have only one input e which drives the valve amplifier; the pilot force is 
represented by an additional state variable. The state equations which govern this plant 
are 



" X P 


-k/c 

A/c 

0 

0 

a k/c 

a 

0 

X P 


0 


P 


4A/Sk/V E C 

-4/3(A 2 + C p c)/v E c 

0 

0 

-4A«^k/v E C 

-4 A crj9/V E 

0 

P 


4^!C 1 JC 2/ /V e 


6 


0 

0 

0 

1 

0 

0 

0 

6 


0 

_d 

6 

_ 

0 

0 

-G/I 

-B/I 

0 

0 


6 

+ 

0 

dt 

x b 


0 

0 

0 

0 

0 

1 

0 

x b 


0 


x b 


0 

c*R b 2 A/j 

0 

0 

0 

0 



I 

0 


F P 


0 

0 

0 

0 

0 

0 

0 

F P 


0 


(19) 


In this model, the flow control valve dynamics are neglected. This approximation was 
found to be reasonable for design purposes and advantageous when implementation of the 
compensation is considered. As equation (19) shows, the pilot input force Fp can be 
described by the differential equation dFpy/dt = 0, whose solution for Fp(0) * 0 is a 
constant. Thus, in the compensation design procedure, Fp is a step function. Note that 
the analog computer is included in the plant description. 

_T- 

The compensation design was accomplished by computing a feedback law e = c x 
and performance index V so that 


V = 



( 20 ) 
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where [qJ is a positive definite matrix, R 2 is a positive scalar, c is a constant vec- 
tor, and x is defined as in equation (19). The matrix [qJ was chosen so that the qua- 
dratic term x^[jq]x became 

x T [Q]x = 100(6 - 9) 2 + Rj(6 - 9) 2 (21) 

Therefore, the performance criterion was the minimum integral of the squared difference 
between desired and actual output displacements and velocities. The term R 2 e2 in 
equation (20) provides a linear limit constraint on the valve amplifier input e. The max- 
imum magnitude of e, and hence maximum valve current and flow rate, can be controlled 
by adjusting the weighting factors Rj and R 2 . 

The feedback coefficient vector c was computed by approximating the integral in 
equation (20) by a sum and applying a discrete dynamic programing algorithm (ref. 4). 

The numerical algorithm resulting from these operations is 

i. - -|r 2 + [H] T [E]._ 1 [ H ]j [Hf [e]^*] (22a) 

[■], - ([•] * ♦ [H]5j) + [Q] ♦ Cj T R2Cj (22b) 

The determination of the optimal feedback c is accomplished by assuming that 

[E] = [q] and then recursively computing [e] and c until the difference Cj - Cj_j is 

arbitrarily small. The weighting factors Rj and R 2 must be chosen so that the actual 
plant output is satisfactory (in terms of comparison with desired output) and consistent 
with the physical limitations of the plant components. The fidelity of the actual output 
ordinarily improves as the constraints on the plant are relaxed. In the system considered 
here, the constraint parameters Rj and R 2 are chosen so that the column motion 6 
is qualitatively close to the desired motion 6 for the current and flow rate limits of the 
flow control valve. 

A diagram of the structure of the compensated system is shown in figure 13. Typi- 
cal predicted responses of the compensated system are shown in figure 14, and the numer- 
ical values chosen in the compensation of the system are given in the appendix. The 
responses were simulated by use of the original eighth-order compliant model or equa- 
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tion (17) with feedback from the six state variables included in equation (19). Omission 
of valve dynamics in the compensation procedure did not affect the validity of the feed- 
back law and, consequently, it is not necessary to be able to measure the fluid flow rate 
q and its derivative q for effective control of the system. Computed values of the 
feedback gains c appear to be reasonable for implementation purposes. All the six 
state variables in equation (19) may be measured and conditioned by conventional trans- 
ducers and associated equipment. Adjustments of the feedback gains for variations in 
analog computer settings are required for only the computer output 6 and its. derivative 
6. These adjustments can be accomplished by manipulating input and output gains in the 
computer. The compensated system model demonstrates fidelity over a wide range of 
computer settings and, most notably, eliminates the instability associated with low values 
of the modeled spring rate. A set of typical system parameters and feedback gains cor- 
responding to the plant model examined herein is given in the appendix. Numerical 
experiments have indicated that this feedback scheme is reasonably insensitive to per- 
turbations in plant parameters. 


CONCLUDING REMARKS 

A linear state space mathematical model and state feedback appear to provide sat- 
isfactory means for compensating a hydraulic control loading system containing com- 
pliance in the actuator control linkage. The resulting compensated system can be phys- 
ically constructed by use of conventional components, and it can be adjusted to simulate a 
wide variety of real aircraft by altering input and output gains at the analog computer 
when aircraft parameters are programed. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., October 7, 1974. 
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APPENDIX 


NUMERICAL VALUES 
Physical Constants of Model 

These are the numerical values used in both the rigid and compliant (plant) models. 

A = 15.22 cm 2 (2.36 in 2 ) 

B = 3434.74 N-cm-sec (304 lb-in-sec) 
c = 28.372 N-sec/cm (16.2016 lb-sec/in.) 

C i = 89.309 cm3/sec-mA (5.45 in3/sec-mA) 

C p = 0.24243 cm5/N-sec (0.0102 in 5 /lb-sec) 

F p = 427 N (96 lb) 

G = 108 466 N-cm/rad (9600 lb-in/rad) 

G max = 325 396 N-cm/rad (28 800 lb-in/rad) 
i r = 15 mA (15 mA) 

I = 677.91 N-cm-sec 2 (60 lb-in-sec 2 ) 

J = 422.56 N-cm-sec 2 (37.4 lb-in-sec 2 ) 
k - 10 695 N/cm (6107.86 lb/in.) 

Kj. - 0.0011240 V/N (0.005 V/lb) 

K 2 = 1.56 mA/V (1.56 mA/V) 

K 3 = 3.0047 V/cm (7.632 V/in.) 
p s = 827.28 N/cm 2 (1200 lb/in 2 ) 
q r = 1638.7 cm 3 /sec (100 in 3 /sec) 

R^ = 25.4 cm (10 in.) 

Rp = 79.6925 cm (31.375 in.) 

V E = 191.73 cm 3 (11.7 in 3 ) 
a = 0.59 (0.59) 

/3 = 172 350 N/cm 2 (250 000 lb/in 2 ) 

< = 0.9 (0.9) 

w n = 377 rad/sec (377 rad/sec) 
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APPENDIX — Concluded 


Compensation of the Plant 

These are the numerical values chosen in the compensation of the plant. 

The following values of [q] and R 2 were used for computing feedback coefficients 
according to equation (20): 



0 

0 

0 

0 

0 

0 
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0 

0 

0 
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0 ( 


0 

0 

100 

0 

-10 

0 ( 

[«]- 
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0 

0 

0 

0 

0 ( 


0 

0 

-10 

0 

1 

0 ( 
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0 

0 

0 

0 

0 < 


0 

0 

0 

0 

0 

0 ( 

This [Q] corresponds 

to 

= 0 

in 

equation (21). 


R 2 = 0.1 


Figures 10, 11, and 14 show typical frequency and time responses associated with 
computer settings for a 2.0 hertz, 0.2 damped response of the control column. The feed- 
back coefficient, vector c computed from equation (20), using the values given for [q] 

and R 2 , is given in S.I. Units as 


c = 


-1.94 V/cm 
-0.0000377 V-cm 2 /N 
+22.02 V/V 
+0.669 V-sec/V 
+0.199 V/cm 
-0.00708 V-sec/cm 
-0.00146 V/N 


and in U.S. Customary Units as 


c = 


-4.93 V/in. 
-0.000026 V-in2/lb 
+22.02 V/V 
+0.669 V-sec/V 


+0.506 V/in. 
-0.018 V-sec/in. 
-0.0065 V/lb 
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(a) Gain. (b) Phase. 

Figure 7.- Open-loop frequency response from experimental data. 








Frequency, Hz 

(a) Gain. 

Figure 10.- Overall system open-loop frequency response, 

B = 343.47 N-cm-sec (30.4 lb-in-sec); G = 10 847 



Frequency, Hz 

(b) Phase. 

I = 67.79 N-cm-sec^ (6 lb-in-sec^) 
N-cm/rad (960 lb-in/rad). 


6 and 0 normalized with respect to steady-state value 



Figure 11.- Time response to step force input. 
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I * 677.91 N-cm-sec* 1 
B = 3434.74 N-cm-sec 
G * 108 466 N-cm/rad 


Time, sec 


(b) Compliant model. 
Figure 11.- Concluded. 
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